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We study time-dependent coupled-cluster theory in the framework of nuclear physics. Based on 
Kvaal's bi- variational formulation of this method [S. Kvaal, |arXiv:1201.5548] , we explicitly demon- 
strate that observables that commute with the Hamiltonian are conserved under time evolution. 
We explore the role of the energy and of the similarity-transformed Hamiltonian under real and 
imaginary time evolution and relate the latter to similarity renormalization group transformations. 
Proof-of-principle computations of 4 He and 16 in small model spaces, and computations of the 
Lipkin model illustrate the capabilities of the method. 
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I. INTRODUCTION 

Coupled-cluster (CC) theory was introduced in nuclear 
physics by Coester and Kiimmel [TJH] more than 50 years 
ago, and Cizek and Paldus [3H5] further developed the 
method for applications in quantum chemistry. For re- 
views of this method we refer the reader to Refs. [BHS]. 
The popularity of coupled-cluster theory is due to the at- 
tractive compromise the method offers between accuracy 
on the one hand and computational cost on the other 
hand. In nuclear structure and reactions, for instance, 
the method has extended the reach of ab initio compu- 
tations from light p shell nuclei [TU] to medium-mass nu- 
clei dung. 

Time-dependent coupled-cluster (TDCC) theory dates 
back more than 30 years. It was proposed by Hoodb- 
hoy and Negele [13l [14], with the aim to describe nu- 
clear collisions, and by Schonhammer and Gunnarson [15] 
to compute spectral functions. Until now, however, 
nuclear collision processes are usually described within 
time-dependent mean-field methods [T6lI20j . Likewise, 
applications of TDCC in quantum chemistry have only 
been sporadic [5TH25]. For small amplitude oscillations, 
time-dependent coupled-cluster theory leads to the com- 
putation of excited states within linear response the- 
ory, and the method is routinely used [2TJ [551 US] m 
this framework. However, the situation is different for 
large-amplitude oscillations. We believe that this lack of 
popularity is mainly due to conceptual problems regard- 
ing the conservation of energy, and keeping the energy 
a real (and not complex) number [25]. These aspects 
have only been clarified very recently in Kvaal's formu- 
lation [30] of TDCC as a bi- variational theory [31], and 
the result is a real energy that is conserved under time 
evolution. This approach opens the way to many excit- 
ing applications. In this work, we present a formal proof 
regarding the conservation of energy within Kvaal's for- 
mulation of TDCC and also show that observables that 



commute with the Hamiltonian are conserved under time 
evolution. We study in detail the role of the similarity- 
transformed Hamiltonian under real and imaginary time 
evolution and illustrate our results in applications to the 
nuclei 4 He and 16 O in small model spaces and to the 
Lipkin model. 

This paper is organized as follows. In Sect. [H] we 
re-derive TDCC, set it in relation to time-independent 
CC, and discuss formal aspects about conserved quan- 
tities. In Sect. |III| we focus on the energy. We explic- 
itly demonstrate that the numerical time-evolution of 
the bi-variational CC theory preserves the energy. We 
also study the real and imaginary time-evolution of the 
similarity-transformed Hamiltonian. In Sect. |TV| we dis- 
cuss applications of the method. First, we show that the 
Fourier analysis of the time-evolved coupled-cluster am- 
plitudes of the nuclei 4 He and 16 yields excitation spec- 
tra. Second, we study the Lipkin model and compare our 
results for observables with exact results. We finish with 
our conclusions. Technical aspects are relegated to the 
Appendix. 



II. COUPLED-CLUSTER THEORY 



In this Section, we briefly re-derive the equations that 
govern TDCC theory and prove that observables that 
commute with the Hamiltonian are conserved under time 
evolution. We also consider TDCC in the limit of linear 
response theory, and the limit of vanishing time depen- 
dence. 
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A. Time-dependent coupled-cluster equations 

We are interested in solving the time-dependent 
Schrodinger equation 



H(t)\*(t)) = ihdtMt)) 



(i) 



Here H(t) is a time-dependent Hamiltonian. In TDCC 
theory, the many-body wave function results from an ex- 
ponential excitation operator acting on a product state 



Here 



i*>=n 



4 {t) |o) 



(2) 



(3) 



P =i 



is an A-body product state that serves as the reference 
state. The creation operator a p creates a fermion in 
the single-particle state p, while an annihilation opera- 
tor a p would remove a fermion from the single-particle 
state p. The creation and annihilation operators are al- 
lowed to depend on time, and they fulfill the usual anti- 
commutation relations at equal times. The time depen- 
dence of the creation and annihilation operators is [14] 



a t ot(t) = x;<p(*)i?(*))4(*)- 



(4) 



Here \q(t)) = <xt(i)|0) and \q(i)) denotes the time deriva- 
tive of \q(t)). To simplify notation, we will suppress the 
explicit time argument and denote a q (t) as a q and \q(t)) 
as \q). In Eq. ^ the cluster operator S(t) is 

S(t) = s (t) + S t (t) + S 2 (t) + S 3 (t) + ...+ S A (t) . (5) 

Here so(t) is a complex phase while the operators 
Si(t), . . . , SU. (t) generate lp-lh, 2p-2h, 3p-3h, and up to 
Ap-Ah excitations within the reference state |$), with 



S n (t) = 
1 

(n!)2 



(6) 



E 



, a'.ai 



...a; 



i± ...i n ,ai ...a„ 



The functions So(t) and s° 1 " i a "(t) are the unknown time- 
dependent excitation amplitudes. In what follows, we 
also suppress the explicit time argument of the np-nh 
cluster operators S n (t), the phase so(t), and the excita- 
tion amplitudes s^ 1 "'""(£). In addition, we will use the 
indices i,j,k,... and a, b, c, ... to label occupied and un- 
occupied states in the reference ([3]), respectively, and the 
indices p,q,r, ... to label any state. 

We are interested in the time evolution of the clus- 
ter operator 5. For computational efficiency, we employ 
the coupled-cluster with singles-and-doubles (CCSD) ap- 
proximation and truncate 5 at the lp-lh excitation level 



5 = s + Si + 5 2 



SO + E S i a l a i 



ijab 



(7) 



To derive the equations of motion for the amplitudes 
so, sf, and s° b we left- multiply the time-dependent 
Schrodinger equation ([TJ) with the operator e~ s , such 
that 



H\$) 



ihe- s d t e s m 



(8) 



Here, we introduced the similarity-transformed Hamilto- 
nian 



H=e' s He s 



(9) 



For the computation of the similarity transformation B = 
e~ s Be s of an operator B we use the Baker- Campbell- 
Hausdorff expansion, 



e- s Be s = B + [B,S\ + ~[[B,S\,S\ 



2!' 

[[[5,5], 5], 5] 
[[[[B, 5], 5], 5], 5] 



(10) 



For a two-body operator, this expansion truncates natu- 
rally at its fourth nested commutator. For B = 3 t , the 
expansion in Eq.(10l reduces to 



~ S d t e s 



dt 



■S + \[S,S\ 



(11) 



for any truncation of 5. Note that the commutator [5, 5] 
is only nonzero if the single-particle states depend explic- 
itly on time [see Eq. Q]. 

We insert the CCSD ansatz into Eq. left- 
project with the Hermitian adjoints of the states |$), 
= a\di\§), and = a\a\ajai\<&)] and then use 

Eq. (11). In what follows, we limit ourselves to time- 
independent basis functions. The equations of motion 
for the 5 amplitudes within a time-dependent basis are 
presented in Appendix [S] We obtain 



ihso 
ihsf 



(3|*T|$>, 
(<t>%\H\<f>) 



(12) 
(13) 
(14) 



The left-hand sides of these equations are matrix ele- 
ments of the similarity-transformed Hamiltonian H = 
e~ He . The corresponding algebraic expressions can 
be worked out diagrammatically and are well known from 
time-independent CCSD The resulting expressions 
are linear in the matrix elements of the Hamiltonian and 
nonlinear in the cluster amplitudes. It is useful to rewrite 



the equations ( 12 ) to ( 14 1 in operator form as 



ihS = H%. 



(15) 



Here H\ denotes the operator that is obtained from the 
similarity-transformed Hamiltonian H when all but the 
matrix elements on the right hand side of Eqs. ( 12p4 | are 
set to zero. When viewing H as a matrix (See Eq. (45) 
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below), the first column of the matrix of Hi is thus iden- 
tical to the first column of H , while all other columns are 
zero. 

For the computation of observables within TDCC we 
follow Rcf. [30 . Coupled-cluster theory employs a non- 
Hermitian Hamiltonian H and must thus be formulated 
in a bi-variational approach, i.e. the left (bra) and right 
(ket) states must be varied independently. In this formu- 
lation, the energy expectation value is 



E = ($|LH|$) . 
Here L is a linear de-excitation operator 

1 \ - 



L 



ijab 



(16) 



(17) 



and Iq, l l a , and 1% are the amplitudes associated with up 
to 2h-2p de-excitations [SJ. From Eq. ( 16 ), it is clear that 
the left many-body state is given by 

: ($|Le~ s . 



(18) 



To obtain the time evolution of the operator L, we insert 
of Eq. (18) into the time-dependent Schrodinger 
equation, right-multiply by e s and obtain 



- ihlo — 0, 

-mi = (*\lh\*?), 
-ihi%, = (*\lh\$$) 



(19) 
(20) 
(21) 

Again, we can write these equations in operator form as 

— ihL = LH - LH\ , (22) 

and here it is understood that the right hand side, when 
written in matrix form of Eq. ( 45 ) , has nonzero matrix 
elements in all but the first row. Note that Iq = 1 for 
proper normalization. 



Let us explicitly check that the energy ( 16 ) is con- 

~~ r s He s and Eq. (fTHh 



served. From the definition H 
we find 



H 



[H,S] 



(23) 



Thus 



E = 



($\LH\<Z>) 

i 

ft' 



:((*\{LH 



(<f>\LH\<f>) 
LH X )H\®) 



= 0. 



(24) 



Here we used Eq. |22j) and H |$> = Hi |3>). 

We can extend these results to any observable B that 
commutes with the Hamiltonian. We have for the expec- 
tation value (B) = (3>| LB\$) 



(B) = ($\LB\$) 



($|LB|$) 



(25) 



Here, we used B — [i/K)\H,B\ as the generalization of 
Eq. Q, and [H,B] = follows from [H, B] = 0. Note 
that this proof addresses a long-standing problem already 
raised by Monkhorst [2~4] . 

For the computation of observables, one employs the 
normal-ordered one-body and two-body density matrices 



c ■ 



(Pq P )n 



(26) 
(27) 



Here, only connected (C) diagrams contribute. Due to 
the linearity of L, these densities have finite expansions. 



B. Time-independent coupled-cluster equations 

Let us re-derive the well-known time-independent 
CCSD equations from TDCC theory by requiring that 
the cluster operators Si and S2, and the single-particle 
basis do not depend on time. Thus we obtain from the 
Eqs. (p| and (fl4l 



m\H\$) = , 

($i-\H |$) = . 



(28) 
(29) 



These are the time-independent CCSD equations. They 
state that the reference state |$) is an eigenstate of the 
similarity-transformed Hamiltonian H in the space of Ip- 
lh and 2p-2h excitations. Thus the corresponding eigen- 
value is the energy 



E Q = ($|J?|$) 



and Eq. ( 12 ) yields 



s = --Et 
n 

up to an irrelevant constant. 



C. Coupled-cluster equations of motion 



(30) 



(31) 



The equation-of-motion CC methods [251 E2] are com- 
monly used to compute excited states within the CC 
method. For completeness, we briefly re-derive the corre- 
sponding equations from the time-dependent formalism. 
Again, we assume that the Hamiltonian and the single- 
particle basis do not explicitly depend on time. For the 
solutions of Eqs. (Tl3|) and (14), we make the ansatz 



*?(*) = tt + rf(t), 
stf(t) = i# + r#(t), 



(32) 



assume that the time-independent amplitude s tf an d tfj 
fulfill the time-independent CCSD equations (28|29), and 



that the remainders rf and rfj 1 are small perturbations, 
R with \R\ -c \T\. In leading order we thus 



i.e. S = T 

have 



H 



He 1 



Hn 



(33) 
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We insert S = T + R into Eq. (15), employ Eqs. (281, sq, and periodically compute the difference between the 



(29), and the expansion (10), and keep only the terms 



first-order in R. This yields the operator equation 

ihR = ~Ht + [Ht, R] • (34) 
Then, assuming, a simple harmonic time-dependence, 



rf(t) 



ab -iuit 



(35) 
(36) 



we find that the amplitudes and frequencies solve the 
eigenvalue problem 



hurt = (<S>« \{Hr,R]m 
?Ujt$ = ($$\[Ht,R]\$) 



(37) 



These are the well known equations of motion for the 
computation of excited states within CCSD. 



It is also interesting to study the left equation ( 22 1 for 



S = T + R, with \R\, \L\ <C T. Keeping terms linear in 
R and L yields in leading order 



- iHL = L(H T - E 



(38) 



since (i?r)i = Eq. Again we assume a harmonic time- 
dependence 

£(*) = iy ut > (39) 

l%(t) = (40) 
and find the well-known left eigenvalue problem 

LH T = {E + uS)L. (41) 



total energy, according to Eq. (16), and its initial value. 
Figure [T] shows this difference in total binding energy as 
a function of time for the 4 He nucleus with a time step 
of width 0.05 fm/c. The change in energy is very small 
throughout the time-evolution and energy is conserved 
for all practical purposes. In practice, the precision to 
which the method conserves energy has a notable de- 
pendence on the time step width: we obtain changes to 
energy on the order of 10~ 10 MeV using a step width 
of 0.05 fm/c and changes to energy of average order 
10~ 2 MeV for a step width of 1 fm/c. 




FIG. 1: The time evolution of the change in total energy of 
4 He from its initial value. 



III. ROLE OF ENERGY 



In this Section, we focus on the role of the energy 
in TDCC theory. First, we probe the conservation of 



energy ( 24 ) numerically for simple models of light nu- 



clei. Next, we study the time-dependent eigenvalues of 
the similarity-transformed Hamiltonian H and find that 
these cannot be viewed as energies. This opens the way 
to employ imaginary-time evolutions of the similarity- 
transformed Hamiltonian for the computation of the 
ground-state energy. In what follows we employ a low- 
momentum two-body interaction produced by a simi- 
larity renormalization group transformation [33 of the 
N 3 LO interaction from chiral effective field theory [33] at 
the cutoff A = 1.9 fm -1 . We use a small model space 
of four oscillator shells. For the time integration of S 
and L according to Eqs. (15) and (22), we employ the 



fourth-order Runge-Kutta method. 



A. Time evolution of the energy functional 



For the time evolution, we iteratively solve Eqs. ( 15 ) 



B. Time evolution of H 

We are interested in the eigenvalues of the similarity- 
transformed Hamiltonian Q. For a time- independent 
Hamiltonian H, the eigenvalues of H do not depend on 
time if the similarity transformation is carried out ex- 
actly. For an A-body system, a complete model space 
consists of up to Ap-Ah excitations of the reference state, 
but this is computationally not feasible in general. In- 
stead, one chooses a model space of up to 2p-2h excita- 
tions within the CCSD approximation. This truncation 
renders the similarity transformation inexact and it is 
thus no surprise that the eigenvalues of H are not con- 
stant in time. 

To illustrate the effects of the truncations, we con- 
sider three distinct cases. (?) In the case "TD-CCSD 
(2p-2h)" we employ S = Si + S2 in the computation 
of H [i.e. we time evolve the S amplitudes according 



and (22). We set Iq = 1 and neglect the complex phase 



to Eqs. (13) and (14)], and compute the lowest eigen- 
value of H within a basis including up to 2p-2h excited 
states. Clearly, this case has truncations in the model 
space for nuclei with A > 2. (ii) In the case "TD-CCS 
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(2p-2h)," we employ S = Si [i.e. we set S2 — 0, solve 
the time-dependent coupled-cluster singles (CCS) equa- 
tion (13)]. We compute the ground-state energy within 
a basis including up to 2p-2h excited states. Here, the 
similarity-transformation is simpler than in case (i), but 
the model space is incomplete for nuclei with A > 2. (Hi) 
In the case "TD-CCS (3p-3/i)," we again solve the time- 
dependent CCS equation (13), but compute the ground- 
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FIG. 2: (Color online) The time evolution of the lowest eigen- 
value of the similarity-transformed Hamiltonian H for 3 H 
computed for three distinct cases (see text for details). 

Figure [3] shows the eigenvalues of the similarity- 
transformed Hamiltonian of 4 He as a function of time for 
the three cases discussed above. As expected, the energy 
is not conserved in any of the cases. The fluctuations of 
the energy are smallest for case (Hi) which exhibits the 
least truncation. 

The results of this subsection explain clearly why 
matrix elements (28] or eigenvalues of the similarity- 
transformed Hamiltonian cannot serve as conserved ener- 
gies. However, these results also open up the possibility 
to employ imaginary-time evolution as a projection tech- 
nique. 



C. Imaginary-time evolution of H 

Imaginary time propagation of a given quantum state 
projects out the lowest-energy state that has a nonzero 



-25.3 



state energy within a basis including up to 3p-3h excited w 
states \&fjk)- Here, no approximation is made for a nu- 
cleus with A = 3. In what follows we will consider these 
distinct cases for nuclei with mass number A = 3 and 
A = 4. 

Figure [2] shows the eigenvalues of the similarity- 
transformed Hamiltonian of 3 H as a function of time for 
the three cases discussed above. As expected, the energy 
is only conserved for case (Hi). Case (i) and (ii) are dif- 
ferent similarity transformations and yield different H . 
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FIG. 3: (Color online) Same as Fig. |] but for 4 He 



overlap with this state. We set t = 
obtain the evolution equations 



ihr with real r and 



Or 



ab 



>He b |$) , 
- s He s 1$) 



(42) 
(43) 



from Eqs. ( 13 1 and ( 14 ). It is useful to consider the imag- 
inary time evolution of the Hamiltonian directly. The 
corresponding evolution equation is 



d T H(r) = [H(r),d T S] . 



(44) 



Eq. (44) shows that d T S generates the similarity trans- 
formation of the Hamiltonian. It is similar to the simi- 
larity renormalization group equations employed in tam- 
ing interactions |35H37j and in computations of finite nu- 
clei [38]. The evolution stops (i.e. it reaches a fixed 



point) once the right hand side of Eq. (44) vanishes. At 



the fixed point the matrix elements of H fulfill the time- 
independent CCSD Eqs. ( 28|29 |. Thus, the imaginary 
time evolution emerges as an alternative way of solving 
the coupled-cluster equations. 

As an example, we consider the imaginary time evo- 
lution for 4 He. For the initial conditions we choose 
sf = = s?„ b . We then propagate sf and sf^ accord- 



ing to Eqs. (42) and (43) while periodically computing 



the ground-state energy as the right eigenvalue of H. 
Figure [I] shows the difference AE = \E(t) — £?ccsd| 
as a function of r, obtained using 400 steps of width 
0.5 fm/c. Here, -Eccsd is the energy that one can obtain 
from solving the CCSD equations ( |28|29[ ) directly. The 
convergence to the ground state energy is exponentially 
fast, and the exponent is related to the energy gap to the 
first excited state. 

Let us also visualize the evolution of the r-dependent 
similarity-transformed Hamiltonian. We present the 
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as a function of imaginary time r. 



Hamiltonian in a block format as 



/ ($|.ff|$) 



H 



l$\H\$i) 



|*#) 



(45) 



\{$$\H\$) (^\H\<P%) <*#|ff|*#)/ 



and plot logarithms of the averages of the magnitudes of 
all elements within a given block. Figure [3] shows contour 
plots of these averages taken at r = fm/c and r — 
200 fm/c, with a step size of St — 0.5 fm/c. We see 
that the Hamiltonian is driven to a form that makes the 
reference state |$) a right eigenstate and the element 
the corresponding eigenvalue. At t — 200 fm/c 
the time-independent CCSD Eqs. ([28) and j29] )) have 
practically been solved. 



IV. APPLICATIONS 

In this Section we study two applications of the TDCC 
method. First, we compute the spectra of excited state 
energies for 4 He and 16 O as the Fourier transforms of 
the time-evolved cluster amplitudes and compare the re- 
sults to those obtained by solving the time-independent 
coupled-cluster equations of motion. Second, we revisit 
the interacting Lipkin system and compare with the re- 
sults by Hoodbhoy and Negele [15] . 



A. Energy spectra from Fourier transforms 

In coupled-cluster theory, the standard approach to 
excited states [231 152"] is the solution of the right eigen- 



value problem ( |37[ ). Alternatively, we can obtain a spec- 
trum of excited-state energies from a Fourier transfor- 
mation of any of the time-evolved CCSD amplitudes sf 



and s° b , utilizing the solution of Eq. (|15|. This approach 



involves first solving the time-independent CCSD equa- 
tions, Eqs. (28 1 and (29). Then, we perturb the result- 



ing cluster amplitudes by a small random contribution of 
the order w 0.001 and solve Eqs. (13) and (14). During 



the iterations, we record the time evolution of a single, 
randomly-selected amplitude s$ or . Results for the 
Fourier transform of the CCSD amplitudes for 4 He and 
16 are shown in Fig. [6] and Fig. [7] respectively. The 
length of the time evolution is 20000 fm/c for 4 He and 
6000 fm/c for 16 0, and the step size is 1 fm/c. Cor- 
responding to the length of the time evolution, the en- 
ergy uncertainty in the resulting Fourier spectra is about 
SE « 0.06 MeV for 4 He and SE sa 0.2 MeV for 16 0. In 
each plot (the excitation energies are measured relative to 
the ground state) the giant peaks at E = correspond to 
the ground-state of each system. Note that several of the 
excited-state energies shown in Figs. [6] and [7] are not asso- 
ciated with physical states of 4 He or "O, respectively but 
arise from the limitation of a small model space, and the 
incomplete separation of the center-of-mass motion |39| . 
In Tables [I] and |n] for 4 He and 16 O, respectively, we com- 
pare the values associated with four selected peaks to the 
values computed using equation-of-motion CCSD. Note 
that all of the peaks in Figs. [6] and [7] are associated with 
energies computed with this method. The agreement is 
good, and deviations are within the uncertainties related 
to the finite time evolution. 

Note that TDCC calculation of excited states via 
Fourier transformation requires much more computa- 
tional cycles than the standard approach. The latter 



requires us to solve the eigenvalue problem (37) once, 
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FIG. 6: (Color online) Fourier transform of a randomly- 
selected cluster amplitude s" or Sy for 4 He. Energies are 
measured relative to the ground-state energy. 



TABLE I: Comparison of select excited-state energies ob- 
tained via time-dependent CCSD (TD-CCSD) and equation- 
of-motion CCSD (EOM-CCSD) for 4 He. Energies are given 
in units of MeV. 





TD-CCSD 


EOM-CCSD 


Ex 


3.04 


3.02 


E 2 


29.45 


29.47 


E 3 


35.52 


35.53 


Ei 


41.10 


41.11 



TABLE II: Comparison of select excited-state energies ob- 
tained via time-dependent CCSD (TD-CCSD) and equation- 
of- motion CCSD (EOM-CCSD) for 16 O. Energies are given 
in units of MeV. 





TD-CCSD 


EOM-CCSD 


Ei 


32.20 


32.29 


E 2 


33.85 


33.91 


E :i 


37.40 


37.60 


E 4 


48.15 


48.12 
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FIG. 7: (Color online) Fourier transform of a randomly- 
selected cluster amplitude s" or s^ 6 for le O. Energies are 
measured relative to the ground-state energy. 



degenerate levels separated by an energy e. In the ref- 
erence state, all N particles occupy the lower level. The 
particles interact via the two-body Hamiltonian 



H = 



pa 



V 

2~ 



E 

pqcr 



a pa a qa a q- 



.(46) 



Here p, q = 1, . . . , N are particle labels, a = ±1 denotes 
the upper and lower level, respectively, and V is the 
strength of the two-body interaction. For weak inter- 
actions V(N — 1) < e, the reference state is JI^Li a j>-|0)> 
and the model is in its "symmetric" phase. At V(N— 1) ~ 
e the model makes a transition to a "deformed" phase, 
and a Hartree-Fock reference state should be chosen for 
V(N-1) > e. 

As in Ref. [T3] we consider two Lipkin models (N = 14 
each with Hamiltonians Hi and H2, respectively) that 
are noninteracting for t < and start to interact at t = 
via 



(47) 



while the TDCC method requires one computation of H 1 
per time step. The real advantage of the time-dependent 
coupled-cluster method lies - of course - in the study of 
time-dependent Hamiltonians. 



B. Excitation energies of interacting Lipkin 
systems 

Hoodbhoy and Negele studied the dynamics of two in- 
teracting, 14-particle Lipkin models [ID] with the TDCC 
method [T3]. In this Subsection, we compare our calcu- 
lations to the earlier results. 

The TV-particle Lipkin model consists of two iV-fold 



The observable of interest is the energy difference 
1 



AE(t) ee -^\H 1+ H 2 \y)(t) 
- (*| Hi + H 2 1*) (0) 



(48) 



of the subsystems. 

In Hoodbhoy and Negele's calculation only 2p-2h ex- 
citations were considered (setting Si =0), and Eq. (14) 
is solved. The calculation of expectation values was per- 
formed in a Hermitian approach based on the density 
matrices 



{p qp ) N = ($| e b {rfq}e b m 



c ' 



(Prspq)N 



$|e H '{pVsr}e H |$) 



(49) 
(50) 




FIG. 8: (Color online) Time-dependent coupled-cluster dou- 
bles (TD-CCD) result for the excitation energy as a func- 
tion of time for two interacting, 14-particle Lipkin systems. 
Hoodbhoy and Negele's results for TD-CCD, time-dependent 
Hartree-Fock (TDHF), and the exact calculation are noted 
with an asterisk and were taken from Ref. 1131. 



Here the subscripts N and the curly brackets indicate 
the normal-ordering of the creation and annihilation op- 
erators. The subscript C indicates that the discon- 
nected diagrams have been factored from the expres- 
sions. Whereas these expressions can easily be expanded 
in terms of products of the cluster amplitudes sf and sf!-, 
the expansions are non-terminating due to the exponen- 
tial nature of the excitation operators. Thus they must 
be truncated at some level, and we choose a truncation 
at the terms quadratic in the s^ b amplitudes. Setting 
e = 1 and V = 0.357 for strong coupling, we computed 
the excitation energy of each system as a function of time 
for 12 x 1CP 23 s, making 90 time steps of size 0.4 fm/c. 
Our result for the excitation energy and that obtained 
by Hoodbhoy and Negele [T3] are shown in Fig. [8 The 
results arc indistinguishable. Figure [8] also shows the ex- 
act results and results from time-dependent Hartree-Fock 
(both taken from Ref. [T3"]L 

Let us finally also compute time-dependent expecta- 
tion values within the coupled-cluster framework and em- 
ploy the evolution of the left dc-excitation operator L 
and the density matrix ( 26 ) of the 14-particle Lipkin sys- 
tem. The parameters of the Lipkin model are e = 1 and 
V = 0.04. Thus, (N - 1)V « 0.52 < 1, and we are in 



the symmetric phase of the two-level Lipkin model. We 
choose the initial conditions S = = L and study the 
time evolution of the operator 



Jz 



1 



aal„a v 



2 / j per "per- 

per 



(51) 



Our result shown in Fig. [^exhibits a good agreement with 
the exact solution. The main period and amplitude of 



FIG. 9: (Color online) Time-dependent CCSD evolution of 
the one-body interaction of the 14-particle Lipkin system. 
The exact result is also shown. 



the oscillation is well reproduced but the coupled-cluster 
calculation misses a fine modulation of the amplitude. 
We also performed computations for weaker interaction 
strengths V and observed that the agreement between 
TDCC theory and the exact solution improves with de- 
creasing values of the interaction strength. 



V. CONCLUSIONS 

We have studied the bi- variational formulation of time- 
dependent coupled-cluster theory and illustrated some of 
its features in applications to simple models of the atomic 
nuclei. Our main results are (i) the explicit demonstra- 
tion that observables that commute with the Hamilto- 
nian are conserved under time evolution, (ii) the explana- 
tion of the role of the similarity-transformed Hamiltonian 
under real and imaginary time evolution, and (iii) the 
computation of energy spectra and time-dependent ob- 
servables and their comparison with exact results within 
simple models. These results open the way for the de- 
scription of time-dependent phenomena such as nuclear 
collisions within more realistic models and beyond time- 
dependent mean-field methods. 
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Appendix A: Equations (12), ( 13 | , and (14| in a 



time-dependent basis 

If the single-particle states are allowed to time-evolve, 
the commutator in Eq. (Ill is nonzero, and the equations 



of motion for the cluster amplitudes become 



(&\H\Q) 



= ih 



(Al) 
(A2) 



k kc 

ih\s1 + (a\i) 

c k 

+ 2E s ^( (c| ^- (fc|6) ) 

kc 

kc 

(*>tj\Hm = ih[i^ + ^(8^(a\c)-8^(b\c)) (A3) 

C 



Ipiab^slsfA^k)* -{k\c 



kc 



Here it is clear that the creation and annihilation oper- 
ators are time-dependent, and the permutation operator 
P is defined as P(pq)f(p,q) = f(p,q) - f(q,p)- 
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